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New analysis on narrow baryon resonance decaying into pK 
^-interactions at 70 GeV/c with SVD-2 setup. 



in 



oo 
o 
o 



x 

<D 

Oh 

(N 
> 

7— I 

m 
cn 

o 

oo 
O 



A. Alccv 3 , E. Ardashev 2 , V. Balandin 3 , S. Basiladze 1 , S. Berezhnev 1 , G. Bogdanova 1 , I. Boguslavsky 3 , V. Bychkov 3 , 
V. Ejov 1 , G. Ermakov 1 , P. Ermolov 1 , N. Furmanec 3 , V. Golovkin 2 , S. Golovnia 2 , S. Gorokhov 2 , I. Gramenitsky 3 , 
N. Grishin 1 , Ya. Grishkevich 1 , D. Karmanov 1 , A. Kholodenko 2 , A. Kiriakov 2 , N. Kouzmine 3 , V. Kozlov 1 , 
E. Kokoulina 3 , V. Kramarenko 1 , A. Kubarovsky 1 , I. Kudryashov 1 , V. Kuzmin 1 , E. Kuznetsov 1 , G. Lanshikov 3 , 
A. Larichev 1 , A. Leflat 1 , M. Levitsky 2 , S. Lyutov 1 , M. Merkin 1 , A. Minaenko 2 , G. Mitrofanov 2 , V. Nikitin 3 , 
S. Orfanitsky 1 , V. Parakhin 2 , V. Petrov 2 , V. Peshekhonov 3 , A. Pleskach 2 , V. Popov 1 , V. Riadovikov 2 , V. Ronjin 2 , 
I. Rufanov 3 , D. Savrina 1 , V. Senko 2 , N. Shalanda 2 , M. Soldatov 2 , L. Tikhonova 1 , T. Topuria 3 , Yu. Tsyupa 2 , 
A. Uzbyakova 1 , M. Vasiliev 2 , A. Vishnevskaya 1 , K. Viriasov 3 , V. Volkov 1 , A. Vorobiev 2 , A. Voronin 1 , V. Yakimchuk 2 , 
A. Yukaev 3 , L. Zakamsky 2 , V. Zapolsky 2 , N. Zhidkov 3 , D. Zotkin 1 , S. Zotkin 1 , and E. Zverev 1 

1 Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, 1/2 Leninskie gory, Moscow, 119992 Russia 

2 Institute for High-Energy Physics, Protvino, Moscow region, 142284, Russia 

3 Joint Institute for Nuclear Research, Dubna, Moscow region, 141980, Russia 

March 31, 2008 

Abstract. The inclusive reaction pA — » pK® + X was studied at IHEP accelerator with 70 GeV/c proton 
beam using SVD-2 detector. Two different samples of K®, statistically independent and belonging to 
different phase space regions, were used in the analyses and a narrow baryon resonance with the mass 
M — 1523±2(stat.) ±3(syst.) MeV/c 2 was observed in both samples of the data. The combined statistical 
significance was estimated to be of 8.0 (392 signal over 1990 background events). Using the part of events 
reconstructed with better accuracy the width of resonance was constrained to F < 14 MeV/c 2 at 95% C.L. 
The xf distribution was found to have a peak at zero with < \xf\ >~ 0.1, that qualitatively agrees to a 
Regge-based model predictions. A new cross section estimate of a ■ BR(0 + — > pK°) — 4.9 ± l.0(stat.) ± 
1.5(syst.) fib/nucleon for xf > was obtained. 



• i-H ■ 1 Introduction. 

X' 

' In the last years the observation of a narrow baryon state 
' named <9 + predicted by Diakonov, Petrov and Polyakovp] 
has been reported by a large number of experiments in 

However, several experiments, mostly at high energies, 
did not confirm the existence of + . The complete list 
of references to positive and negative results with dis- 
cussion can be found in P^lfTMTl] . The situation became 
more intriguing when CLAS collaboration reported neg- 
ative results on + photoproduction off proton and deu- 
tron with high statistics [15]. Meanwhile LEPS collabora- 
tion reported on another evidence of 6> + -baryon in the 
reaction -> 6>+y1*(1520)[16]. The STAR collaboration 
observed the double charged exotic baryon in the pK + 
decay channel [17]. Therefore, pentaquark existence is still 
under the question and new experiments are needed to 
confirm or refute it. 

The SVD-2 collaboration has reported earlier the ob- 
servation of narrow baryon resonance in the pi^-system 
with the mass of M = 1526 ± 3(stat.) ± 3(syst.) MeV/c 2 
and r < 24 MeV/c 2 \§\. In that analysis the K° s -mesons 



decayed inside vertex detector were used (decay region 2 - 
35 mm, zero stands for the center of the 1st target plane). 
In this paper we present our new results of the study of 
the same reaction: 

pN -> 0+ + X, 0+ -> P K°, K Q S -> 7T+7T-. 

We have used two independent data samples, selected by 
the point of decay: inside and outside the vertex detec- 
tor ( decay regions 0.2 - 35 and 35 - 600 mm, respectively). 
First results of this work were published earlier f [TSlll9j ). 



2 SVD-2 experimental setup. 

A detailed description of SVD-2 detector and the trigger 
system can be found elsewhere [9"I l2T)ir2"T] ■ For these studies, 
we used the following setup components (Fig[T]): 

1. the high-precision microstrip vertex detector (MSVD) 
with five active (Si) and two passive (C,Pb) target 
planes (AT); 

2. large aperture magnetic spectrometer (LAMS) : multi- 
wire proportional chambers (MWPC) with Y(vertical), 
U{Y - 10.5°) and V(Y + 10.5°) wire directions; 
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Fig. 1. SVD-2 layout (described in the text). 

3. the multicell threshold Cherenkov counter (TCC), with 
the ir/K/p thresholds of 3/11/20 GeV/c. 

Data taking was performed in the 70 GeV/c proton 
beam of IHEP accelerator with an intensity of sw (5 -j- 6) • 
10 5 protons/cycle. The total statistics of 5 • 10 7 inelastic 
events was obtained. The sensivity of this experiment for , 
inelastic piV-interactions, taking into account the trigger 
efficiency, was 2.3 nb~ x . 



3 Sample I: decaying inside the vertex 
detector (decay region 0.2 - 35 mm). 

After the first SVD-2 pentaquark publication[9j (January, 
2004) new improved algorithms of the tracks reconstruc- 
tion had been developed. It increased 3-4 times the avail- 
able number of K® decays inside the vertex detector and 
essentialy improved the experimental K® mass resolution. 
The selected events were having a well defined secondary 
vertex in the region of 0.2 - 35 mm along the beam direc- 
tion from the beginning of the active target (correspond- 
ing to the sensitive area of the vertex detector). Secondary 
tracks were traced to the magnetic spectrometer to obtain 
their momenta. About 3T0 4 kaon candidates in the events 
with the number of primary charged tracks N c h < 8 were 
selected, to reduce the combinatorial background in fur- 
ther analysis. Primary particles were traced to the spec- 
trometer and two-body invariant mass distributions were 
built using previously selected neutral particles in asso- 
ciation with primary vertex tracks. K®ir + and Air + in- 
variant mass spectra showed that the masses and widths 
of K* + (892)-meson and I7 + (1385)-baryon are consistent 
with PDG tables [22]. 

The pK® invariant mass spectrum is shown in Fig [2] 
An excess of s=205 signal over b=1050 background events 
is observed in the area of interest, with the significance of 
s/V^+b = 5.8 . 
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Fig. 2. Sample I: The (pK®) invariant mass spectrum for the 
events with N charged < 8. The dashed line stands for the peak 
extracted from the fit. 
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Fig. 3. Sample I: The {pK®) invariant mass spectrum for 
decaying inside the vertex detector with additional quality cuts 
explained in text. 



To estimate the natural width of the observed peak 
we selected events with the best reconstruction accuracy, 
applying the track quality cuts: 

- D <2a, 

where D is the distance of the closest approach for V° 
tracks and a is a calculated error, 

- N hits > 12, 

where Nhits is the number of hits for the track in the 
spectrometer planes. 
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Fig. 4. Sample II: a). The (tt + -k~) invariant mass spectrum. 
A window corresponding to ±3u is shown by the vertical lines, 
b) The (pn~) invariant mass spectrum. 



With that, we obtained the mass resolutions 3.98 MeV/c 2 
for K° and 1.46 MeV/c 2 for A particles. 

The background was reduced further with the selection 
of events with 3 GeV/c < P pro ton < 10 GeV/c and N c h > 
3; the latter cut may correspond to the greater inelastic- 
ity of the events with the + creation [23]. The result is 
shown in Fig[3j Here the curve stands for the normalized 
mixed-event background built of kaon and proton taken 
from different events with the same kinematic require- 
ments. The distribution was fitted by the sum of Gaus- 
sian function and second-order polynomial background. 
The Gaussian sigma of 4.1 ± 1.0 MeV/c 2 is consistent 
with our estimate of the highest possible experimental res- 
olution of the SVD-2 setup. Taking into account the ex- 
perimental SVD-2 resolution (calculated using well-known 
resonances) it was estimated that the intrinsic width of the 
(jpK° s ) -resonance is r < 14 MeV/c 2 at 95% C.L. No signal 
was observed in the combinations of K® with a tracks of 
negative charge. 

4 Sample II: "Distant" decaying outside 
the vertex detector (decay region 

35 — 600 mm). 

To create a Sample II, -fT^-mesons and A-hyperons were 
reconstructed by the following procedure. The events with 
n c h < 8 in the primary vertex were reconstructed using 
the algorithm described in [§] and the spectrometer hits 
belonging to the reconstructed tracks were removed. The 
remaining spectrometer hits were used to search for tracks 
originating from secondary vertices in the region before 
the first spectrometer plane (decay region 35 — 600 mm). 
The initial track candidate had to have at least 3 hits in 
the Y-planes which could be approximated by a polyno- 
mial function originating from the area of interest. For 
these candidates the hits in the U and V planes of the 
spectrometer were searched for and global track param- 
eters were defined using magnetic field map. Tracks of 
opposite signs were combined to test on having a common 
secondary vertex. Based on the intrinsic tracking resolu- 
tion of the spectrometer the minimum distance between 
two tracks was required to be less than 1 mm in the hori- 
zontal and 5 mm in the vertical directions. The efficiency 
of this method was estimated as 50%, it's lower than that 
for the reconstruction based on vertex detector informa- 
tion. 

The resulting invariant masses of (7r + 7r~) and (p7r~) 
combinations are shown in Fig[4^, and |4|d. The standard 
deviations in the mass distributions are 6.1 MeV/c 2 and 
2.4 MeV/c 2 for K® and A masses respectively. The ob- 
tained tracks collection was used as an opportunity to 
check for the production of well-known resonances. For 
</>°(1020) -> K+K- and A* (1520) -► K~p (FigE} a TCC 
particle identification was used. The total number of de- 
tected /1*(1520) — ► K~p decays was about 2 ■ 10 4 . Using 
PDG data, we estimated our experimental mass resolu- 
tions as cr y i»( 152 o) = 2.4 ± 0.9 MeV/c 2 and cr^o^ W2 a) = 
1.6 ±0.2 MeV/c 2 . 



The production of the resonances decaying into strange 
neutral particles was also studied. The (ir + K®) and (Air + ) 
invariant mass spectra are shown in Fig|6] The clear peaks 
from 7^* (892) and .27 + (1385)-resonances are seen. In total, 
the masses and widths of well-known resonances were all 
in a good consistency with their PDG values (|22j). For 
a 17+ (1385), when applying a pa < 6 GeV/c cut, we ob- 
served a structure near 1480 MeV/c 2 (Fig(7j) . This peak 
may correspond to the 17(1480), marked as one star reso- 
nance in the PDG review f[2"2"]). 

To search for the <9 + -particle, the events with multi- 
plicity of six or less charged tracks in the primary vertex 
were selected to minimize the combinatorial background. 

About 52000 if °-mesons in the mass window of ±20 M eV/c 2 
were found in the selected events. To eliminate contami- 
nation from A decays, candidates producing the effective 
mass of less than 1.12 GeV/c 2 for mass hypothesis 

were rejected. The resulting invariant mass distribution is 
shown in Fig0^. 

Proton candidates were selected as positively charged 
tracks with a number of spectrometer hits more than 12 
and momentum 8 GeV/c < P p < 15 GeV/c. Pions of such 
energies should leave a hit in the Threshold Cherenkov 
Counter, so the absence of hits in TCC was also required. 
Moreover, for the events with a TCC hit of a positive par- 
ticle taken as a proton, no significant peaks were found 
in the (pK®) invariant mass spectra. We estimated that 
the 7r + -background averages to no more than 10% under 
selection criteria used. 
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Fig. 5. Sample II: a) The (K + K ) invariant mass spectrum, 
b) The (K~p) invariant mass spectrum. 
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Fig. 6. Sample II: a) The (ir + Kg) invariant mass spectrum, 
b) {A-k + ) invariant mass spectrum. 



Effective mass of the pK® system is plotted in FigJS] A 
peak is seen at the mass M = 1523.6 ± 3.1 MeV/c 2 with 
a a = 12.9 ± 2.5 MeV/c 2 . The distribution was fitted by 
a sum of a Gaussian function and second-order polyno- 
mial background. There are s=187 events in the peak over 
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Fig. 7. The (An + ) invariant mass spectrum with the addi- 
tional pa < 6 GeV/c cut. 
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Fig. 8. Sample II: The (pK^) invariant mass spectrum with 
the cuts explained in text. The dotted line in the bottom stands 
for the peak extracted from the fit. 



b=940 background events. The statistical significance for 
the peak is s/y/s + b = 5.6 

Two different models were tried to describe the back- 
ground. The first was taken from RQMD Monte Carlo 
(FigJS^,) [53] and the second from the mixed event ap- 
proach (FigJSJa). Both of them fitted plausibly the exper- 
imental data. 
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Fig. 9. Sample II: The (pK°) invariant mass spectrum with 
different background descriptions represented by fitted dashed 
histogram: a) mixed-event model background; b) RQMD 
Monte Carlo background. 
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Fig. 10. The (pK®) invariant mass spectrum for positively 
charged particles considered as a proton obtained from RQMD 
simulation. Dashed histogram represents the (pKg) invariant 
mass spectrum with identified protons. Inset shows the (nK®) 
invariant mass spectrum. 



It was verified that observed p-fT^-resonance can not 
be a reflection from other (for example if** (892) or A ) 
resonances. The mechanism of producing spurious peak 
around 1.54 MeV/c 2 involving K® and A decays overlap 
as suggested by some authors [26,27 was also checked and 
no ghost tracks were found. 

5 Monte-Carlo simulation of pA-interactions. 

The invariant mass spectra of different particles and 
in pA collisions have been simulated by means of RQMD 
2.3 event generator [23]. 

The RQMD generator produces spectra including short 
lived resonances. We used these spectra as input data for 
GEANT 3.21 package [23] to achieve a set of final state par- 
ticles (K° and others) which are detected by apparatus. 
This set of particles was used to implement experimental 
cuts and to produce invariant mass spectra. The (uK®) 
invariant mass spectrum clearly shows peak of yl*(1520) 
with 30% enhancement over combinatorial background. In 
contrast, the mass spectrum of the pK® system is smooth 
and contains no fluctuations exceeding statistical errors 
on the background(Fig flOj) . 

6 The x^-distributions and cross sections. 

We performed two searches for + particle using SVD 
'2002 data. Their independence is based primarily on us- 



ing V candidate vertices in different regions. In both 
cases, numerous checks were made on possible peak imi- 
tations. For the kinematic reflections (K*(892) + — > K°p, 
£(1385) + — > K°p, etc.) and mechanisms with ghost tracks 
described in [2 6 [127] null results were obtained. A smooth 
shape of mixed-event background also excludes the possi- 
bility of generating the sharp peak due to the applied cuts 
and the detector acceptance. A narrow shape of the peak 
by itself makes it very difficult to be generated by any 
kind of reflections or kinematic constraints on the data. 

It is impossible to determine the strangeness of this 
state in such an inclusive reaction, however we did not 
observe any narrow structure in (An + ) invariant mass 
spectrum in 1500 ~ 1550 MeV/c 2 mass region (Figj6p), 
expected for a particle with the negative strangeness. 

Our search for (9 + -particle is an inclusive experiment 
with a significant background contribution. We have made 
an attempt to apply a subtraction method to investigate 
creation region in terms of xp. An effective mass distribu- 
tion around the peak was fitted with a sum of background 
(B) and Gaussian functions, the latter having a mean 
of Too and a standard deviation of a. Background was 
taken as a product of a threshold function and a second 
order polynomial P2, B(m) = P2(m)(l - e"*^™-^)), m 
standing for the effective mass. All the fit parameters were 
given reasonable seeds but no boundaries to prove a fit sta- 
bility. The peak region was defined as too— 2a < m < to,o + 
2a. A number of effective background events under the 
peak, N(B pea k), was evaluated by integrating background 
function over the peak region, x f distributions were taken 
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separately for a peak and an out-of-peak regions ( "wings" ); 
the latter ranged from 1.4 to 1.7 GeV/c 2 with a peak 
region cut out, and a result was scaled to a N(B pea k). 
Assuming that the background characteristics were uni- 
form, we subtracted "wings" distribution from the peak 
one. These operations were performed separately over the 
data from both samples. 

Acceptance corrections were evaluated using a GEANT 
Monte-Carlo [25] for the + registration only. We took into 
account a presence of K° decay in the selected region, 
geometry acceptances, instrumental efficiencies and cuts 
used in the analyses. We did not consider any other tracks, 
accompanying the theta creation. Theta particles were 
generated with a flat distribution over the —1 < xf < 1 
region with a < Pt >= 0.52GeV/c, estimated from a Sam- 
ple I data and typical for the heavy particle creation. A 
vertex of the decay was generated randomly over the mea- 
sured beam spot in X and Y, and Z was taken in the cen- 
ters of the targets. Reconstruction procedure was applied 
to the simulated detector hits, taking into account detec- 
tor efficiencies. The overall result is shown in fig. QT] The 
lower absolute value and a fall at xf — ► 1 for the Sample 
I are dominated by the limited region allowed for the K® 
decay, and the narrow shape for the Sample II reflects the 
proton momentum cut used. 




Fig. 12. xf distribution for the Sample I. Crosses are ex- 
perimental data, curves stand for a quark fragmentation(solid 
line) and a sum of quark and diquark fragmentation(dotted 
line) in a Regge-based model [28]. Y-axis: acceptance-corrected 
cross section over 0.1 in xf- All the data are normalized to the 
4.9 /ib/nucleon for < xf < 0.3 region 
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Fig. 11. Efficiencies of the + registration in xf- Data for 
different targets are weighted using the total numbers of K° 
events. The result for a Sample I is scaled by 10. 



Fig. 13. The same as in fig |12l for the Sample II. 



certain difference in a total number of events. The latter 
makes a contribution to the cross section error calcula- 
tions. 
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The results are shown in figsfT^] and Q2] We plot also 
normalized curves of the predictions made in a Regge- 
based model[28]. In this model, an xp distribution comes 
from the sum of quark fragmentation (bell-shaped around 
Xf = 0) and diquark one (seagull-shaped), each one with 
its proper weight. In [35] the weights were taken as 1:10, as 
coming from the analysis of non-exotic baryons creation. 
Our data may indicate a favoring to the quark fragmen- 
tation part of the model. 

The Sample II has a specific narrow acceptance in if 
due to the proton momentum restrictions. We projected 
the x ^-result of Sample I to a Sample II to check the con- 
sistence of our observations. It was made by applying an 
inverse acceptance correction for Sample II to the result 
of Sample I. It is shown in hg[l4] together with a raw (no 
acceptance correction) data for a Sample II. We found 
a plausible agreement of the distribution shapes and a 



Fig. 14. A comparison of the Sample I (dotted histogram) to 
a Sample II (crosses) results. 

We have made an estimate of A-dependence in the 
+ creation. With SVD-2 vertex detector, Z-coordinate of 
primary vertex was reconstructed to a precision of about 
70 /im, so the events belonging to different target planes 
were easily separated. For each material, a ratio was taken 
of the number of <9 + events to the number of all the 
events. Weighting it to a C 12 value (table [T]), we found 
no differences in the material dependence. Assuming that 
the observed K° events exhibit a /l on -dependence, we 
estimate that + creation follows the same law (within 
errors). 

The cross section of + creation in pN interactions 
for xf > was evaluated from Sample I data (as having 
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Table 1. A-dependence. The fractions of (9 + -events in the K' 
selection for different materials, weighted to a C 12 result. 



Target 


Relative cross section 


c m 


1.0 ±0.22 




0.87 ±0.29 




1.09 ±0.32 



a better acceptance at xp — ► 0): 

ae pN ■ Br(€> + ^pK ) = — — — 

0.34 Jy pN ■ e© 

= 4.9 ± l.0(stat.) ± l.b(syst.) fib/nucleon, 

where <j p n = 31.14 fib is the inelastic proton-proton cross 
section^, 0.34 is the PDG[22] value for the Br(K° -> 
7r + 7r~), Nq and NpN are the numbers of events with a & + 
particle and all the events respectively and e© is the effi- 
ciency of & registration. A luminousity correction factor 
ficorr = 0.84 was evaluated from the following considera- 
tions. A trigger used for the physics events was developed 
initially as a Level I trigger, designed to be sensitive to all 
the inelastic pA-events in the target planes. A compari- 
son to the "beam" trigger events showed that in our sam- 
ples a registration of the primary vertices with a charged 
particles multiplicity N c h < 3 is suppressed. We estimate 
the fraction of total inelastic events accepted by trigger 
as 0.63. Next, in a Sample I we cut out the events with 
N c h > 8 accounting for 0.4 of all the registered events. As- 
suming that + -creation is 1) suppressed in the low N c h 
region and 2) does not depend on the multiplicity in the 
high N c h region, we derived a combined scaling coefficient 
for the <9 + cross section as e cor r — 0.84. As the better es- 
timate of the cross section dependence on the multiplicity 
was not available, the resulting systematic error is rather 
large. With that correction, our integrated luminousity 
was 2.3 • 10 33 cm -2 . It was implicitly assumed that ctq 
in pA interactions also scales as A . This cross section 
value agrees within errors to the result published in [19 , 
but differs from our earliest estimate [9]. The main reason 
is an unexpected form of xp distribution, assumed flat in 
our first publications. Taking into account our situation of 
using inclusive data, and undiscovered yet mechanism of 
theta particle creation, we believe that the cross section 
value is subject to future investigations. 

As positive as negative results of the theta particle 
searches were very often presented as the ratio of theta 
to 4(1520) cross sections. We estimate the latter value as 
147 ±35 fib/nucleon, that agrees with the result obtained 
by EXCHARM collaboration [30] and gives the ratio of 
0.07 assuming that the <9 + — > pK° and + — > nK + de- 
cays have equal probabilities. However these two particles 
may have quite different creation mechanisms. For exam- 
ple, extracting xf for .4(1520) by the same method as de- 
scribed above resulted in much more forward-oriented dis- 
tribution ffig [T5|) . than for a theta particle. It weakens the 
reasons of making this comparison. We may suppose that 
a best way to present such a ratio would be to accompany 
it with corresponding acceptances and xp distributions. 
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Fig. 15. An x F distribution for a /1(1520). 



7 Summary and Conclusions. 



The inclusive reaction pA — » pK^+X was studied at IHEP 
accelerator with 70 GeV/c proton beam using SVD-2 de- 
tector. Two different samples of K® , statistically indepen- 
dent and belonging to different phase space regions, were 
used in the analyses and a narrow baryon resonance with 
the mass M = 1523 ± 2{stat.) ± 3(syst.) MeV/c 2 was ob- 
served in both samples of the data. We observed the total 
of s=392 signal events over b=1990 of background. Using 
additional track quality cuts we obtained r < 14 MeV/ c 2 
on 95% C.L. The statistical significance of the peak can 
be estimated to be (for different estimators): 
s/V5 = 8.7, s/y/s + b = 8.0, s/V^+2b = 5.9. 
The new analysis and larger statistics confirm our previous 
result [9]. The resonance observed is not a statistical fluc- 
tuation neither induced by background processes. We plan 
to continue our investigations of the creation mechanisms 
of pK® resonance, momentum and angular dependencies 
and cross sections. The cross-section estimate for xf > 
(cr • BR{0 + — > pK°) = 4.9 fib) was obtained, xp distribu- 
tion was found to peak at zero with < \xf\ >~ 0.1, that 
agrees qualitatively to a Regge-based model[28]. While in 
agreement with evidences of + observation, there is no 
direct contradiction to null results in hadron-hadron fixed 
target experiments (e.g. |31j). mainly due to different ac- 
ceptances at xf ~ 0. 
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